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We study the effects of imperfections on the fidelity of the Toffoli gate recently realized in a 
circuit QED setup. The noise is introduced to the interqubits interactions hence the couplings no 
longer remain constant but belong to the interval [1 — 5, 1 + 5] J, where J is the average coupling 
constant and parameter 8 quantifies the noise level. We consider both static and dynamical im- 
perfections. We also obtain a more robust gate by solving the quantum optimization problem with 
a modified objective functional including an interval of coupling values rather than a single ideal 
value. 
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I. INTRODUCTION 

Physical implementations of quantum information pro- 
cessing are always subjected to various imperfections that 
decrease the performance of the process. Dynamical im- 
perfections as a result of system-environment coupling 
produce decoherence in the system destroying the bene- 
fits of using quantum information. Static imperfections 
on the other hand do not introduce decoherence to the 
system yet lead to error as well. 

A two dimensional lattice of qubits with nearest- 
neighbor interqubit couplings has already been consid- 
ered as a standard generic quantum computer model to 
incorporate imperfections [HE]- This model shows that, 
for a system affected by static imperfections, quantum 
chaos sets in above a critical interqubit strength and an- 
nihilates the quantum computer performance [2] . Conse- 
quently, the entanglement dynamics also exhibit a transi- 
tion from integrability to quantum chaos [3J 0]. However, 
the dynamics remain almost unaffected for disorder less 
than 10% [6] . The same model was used to study dy- 
namical imperfections in quantum computers in Ref. [7J, 
where a characteristic frequency is associated with the 
noise specifying the rate at which the noise changes. It 
was shown that for low frequencies the imperfections can 
be considered static, and for sufficiently high frequencies 
the effect of noise completely disappears. 

Implementations of quantum computers specifically 
require high fidelity quantum gates. There are many 
measures for the robustness of a quantum gate against 
noise [Sj . Particularly, one- and two-qubit quantum gates 
that are universal for quantum computing have already 
been analyzed under the influence of noise [9lH2j. How- 
ever, the implementation of multi-qubit gates using stan- 
dard decompositions in terms of a universal gate set may 
not be efficient because the implementation time may 
exceed the decoherence time. It is interesting to analyze 
efficient ways of implementing multi-qubit gates directly. 

The Toffoli gate, a three-qubit gate with central role in 
quantum information processing, has been recently im- 
plemented with less resources and moderate fidelities [i"3l - 
[13] , However, new proposals have been suggested to re- 



alize the Toffoli gate with fidelities above 99% [16] [T7]. 
Specially in Ref. [T5] , a quantum control scheme has been 
proposed to realize the Toffoli gate in circuit quantum 
electrodynamics (circuit QED) set up. The gate time is 
about 140 ns that is fast enough because the decoher- 
ence time T 2 in such systems is 10 to 20 [is. The usage of 
quantum control schemes was analyzed in Refs. [TEl - f2"5] 
and circuit QED in Refs. [21H2T] . 

In this paper, we study the effect of imperfections on 
the Toffoli gate implementation. The Hamiltonian is bi- 
linear with XY-type Heisenberg chain for the system and 
Zeeman-like for the control part. Both parts are sub- 
jected to noise. The effect of noise on the control fields 
has already been considered [T51 [32 US] showing that the 
gate is more robust when the single control pulse dura- 
tion is reduced. Actually, for a fixed gate time, increasing 
the noise has less effect on the average fidelity of the gate 
with higher number of control pulses. However, we in- 
vestigate the effect of both static and dynamical noise on 
the interqubit couplings. We also use a control method 
to find a new set of control fields improving the fidelity 
affected by noise. The results are based on numerical 
simulations. 

The paper is organized as follows: Sec.[H]is an overview 
of the quantum control method proposed in Ref. j!6j . 
which realizes the Toffoli gate in circuit QED. In Sec. |III[ 
we describe the noise model and analyze the effect of the 



dynamical and static noise on the gate fidelity. In Sec. IV 



we obtain a more robust Toffoli gate by modifying the ob- 
jective functional. Finally, the summary and discussion 
are presented in Sec. [Vj 



II. IMPLEMENTATION OF THE TOFFOLI 
GATE 

Transmon qubits can be coupled together within 
circuit QED through a transmission line resonator. 
They can be isolated from electromagnetic environment 
and controlled through resonant microwave drives [28) . 
Quantum control techniques can then be applied in this 
setup to implement various quantum information pro- 
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cessing tasks [T5| . 

The interaction Hamiltonian of an array of three trans- 
mon qubits coupled to a superconducting transmission- 
line resonator can be effectively described by an XY-type 
(flip-flop) Hamiltonian PH H5] 

H = Jij {(Tix<rj x + ViyCTjy) , i,j = 1, 2, 3 (1) 

where Jij are coupling constants and Oi X and CTiy are 
Pauli X and Y matrices, respectively. This system is 
manipulated by a Zeeman-like Hamiltonian 



Hc(t) = 5>< <} + 4° (2) 

i=l 

with control fields Ux (t) and u y (t) affecting the qubits 
in x and y directions. Therefore, the system dynamics is 
governed by the total Hamiltonian 



interval lead to the time evolution operators U x 



H{t)=H + H c (t). 



(3) 



These control fields can be implemented through wave 
generators. However, to keep the transmon qubits well- 
defined two-level systems, the fields cannot be arbitrarily 
large. The norm 



max 

i,t 



aw 



i4 ] 



(4) 



is then restricted to be smaller than some threshold value. 

The controllability issues of the system can be ver- 
ified by considering the Lie algebra generated by 

y } which is actually the 
Lie algebra su(8). This algebra provides the operator 
controllability of the system [TB] . 

Finding the control fields that implement the Toffoli 
gate is a numerical optimization problem. A given gate 
time t g is divided into N t (even number) intervals of the 
same time length T. Then, the set of fields is supposed to 
be constant in each interval and are acting on the related 
qubits in alternation in the x and y directions. 

In the first interval < t < T, three x-control pulses 

(i) 

with constant amplitudes u x \ are applied to all three 
qubits. The dynamics is then governed by the Hamilto- 



H x 
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i=l 



In the second interval T < t < 2T, the y-control pulses 

(i) 

with amplitudes u y \ are applied leading to the Hamilto- 
nian 



exp(—iH Xin T) and U y 



exp(—iH y ^ n T) respectively 



in corresponding intervals. The product U y , n U x ^ n , for 
n = 1, . . . ,N t /2 in the reverse order, is the time evolu- 
tion operator U(t = t g ). 

The values of the 37V f control fields are obtained 
through maximizing the fidelity 



(5) 



where u is the concatenation of all control pulses and 
Uto« is the Toffoli gate, given by 
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C^Toff = 



(6) 



in the computational basis. The local optimized control 
fields are obtained through numerical maximization of 
the fidelity with an initial guess for the fields. The min- 
imum gate time for a given fidelity can also be obtained 
by starting with a sufficiently large t g and then gradually 
decreasing its value. 

The piecewise-constant control pulses obtained here 
can be filtered through a low-pass filter to be possible 
to be generated by an actual wave generator. Using a 
product formula approach (see Ref. [55]), it is possible to 
obtain the fidelities for the filtered control fields. 

Assuming w max < 130 MHz, J = 30 MHz and 
J\2 = J23 = 6J13 = J, the Toffoli gate can be real- 
ized with a fidelity larger than 99% in 4.I8J -1 = 140 ns 
with Nt = 20 control pulses. In the following sections we 
focus on this specific realization. 



III. THE TOFFOLI GATE IN SYSTEMS WITH 
IMPERFECTIONS 

In this section we investigate how the gate fidelity is 
affected by the noise in the system. Therefore the sys- 
tem Hamiltonian Eq. (nl undergoes noise while the con- 
trol Hamiltonian Eq. (f2T) is maintained perfect (without 
noise). The scenario is to apply a set of control fields 
optimized for the perfect system to an imperfect system. 
In other words, we analyze the sensitivity of the control 
fields to the noise in the main system. 
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H Vt i = H Q + ^2u y l \a lv . 

i=l 



This process is repeated N t /2 times to complete all N t 
intervals. These time-independent Hamiltonians for each 



A. Noise Model 

We are specially interested in the case where the inter- 
actions between the qubits are noisy. Each Jy in Eq. |lj) 
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is considered to be independently coupled to a stochastic 
variable described by 



Jij{t) = JijO- + 6ij(t)), 



(7) 



where are average coupling strengths and ey (t) are in- 
dependent Gaussian random variables. This noise model 
was used in many physical contexts and can be seen as 
an effective behavior when we average out the effect of 
the environment over the main system |30j . 

We suppose that all €ij(t) have the same mean and 
variance and change their values in time simultaneously 
with a fixed frequency f c . Then, they remain fixed during 
the time interval r c = l// c : 



e(t) 



(fc - 1)t c <t<kr c 



(8) 



where we have neglected the subindices because all com- 
ponents have the same behavior. They are indepen- 
dent and identically distributed (iid) Gaussian random 
variables in different time intervals with expectation 
E[e( fc )] = and variance E[(e^ fc ^) 2 ] = a 2 < oo, where 
fc runs over integer values. 

The autocorrelation function for e(t) is given by 



(e(t)e(t + r)) 



for t in [— t c ,t c ] and zero otherwise. The power spectral 
density (Fourier transform of the autocorrelation func- 
tion), which displays the essence of the noise, is 



O T c . 2 /WT t 

, sine — — 

2ir V 2 



(9) 



In Sec. IIIIBI we describe how we simulate the effect of 
noise on the fidelity for values of r c such that r c < t g . For 
r c > t g we are only interested in the limit r c —¥ oo. In 
this case the power spectral density function approaches 
to the delta function on zero, which corresponds to a 
fixed noise in the system. It means that the value of the 
random variable remains fixed for ever. Physically, this 
situation is associated with inaccuracies in the system 



parameters. This situation is analyzed in Sec. IIIC 



B. Dynamical Imperfections 

We start by obtaining a set of optimal control fields 
implementing the Toffoli gate in the perfect system. The 
fields can be found by maximizing the fidelity given by 
Eq. J5j) using the parameter values given at the end of 
Sec.[n[ To be close enough to the global optimal solution, 
we carry the maximization process over 200 random ini- 
tial guesses and then select the set of fields with largest 
fidelity 

In our noisy system, the couplings (see Eq. Q) 
remain no longer fixed and evolve according to Eq. (ITl. 
There are three different coupling constants and, in prin- 
ciple, the noise affect each of them independently. We 



are imposing that Jn — J23 = 6J13 = J and we are as- 
suming that the random variables €ij(t) associated with 
those couplings have the same mean, variance and char- 
acteristic frequency. Under those assumptions, we will 
show later on in this section that introducing the noise 
only in J is essentially equivalent to introduce the noise 
in the three coupling constants. This explains why we 
can drop the subindices of eij(t) and call the random 
variable simply by e(t). 

Let e(t) change its value with a fixed frequency f c = 
l/r c . We generate t g f c random numbers independently 
according to a Gaussian distribution with a zero mean 
and standard deviation a. Accordingly, we have a realiza- 
tion of e(t) and therefore J(t) for the whole interval [0, t g ]. 
It is now possible to calculate the time evolution operator 
using the total Hamiltonian Ho(J(t)) + H c (t) with those 
values of J(t) above realized. The corresponding fidelity 
for that realization can then be calculated using Eq. {§}. 

The next step is to calculate the average fidelity, which 
is obtained by repeating the above process with the same 
f c and a for a large number of realizations of e(t) within 
the corresponding time interval and summing over all 
fidelities and dividing by the total number of realizations. 
Finally, using this method systematically, we obtain the 
fidelity as a function of a with a fixed f c . 

Fig. [T] shows the average fidelity as a function of a for 
t g fc = 200, 100, 40, 20, 10, 5, 2, 1 from top to bottom, re- 
spectively. The number of realizations per each fixed a 
for the first four items is 10,000 and for the last four items 
is 100,000. The average fidelity decreases for all frequen- 
cies when the standard deviation a increases, confirming 
our intuition about the effect of the decoherence over 
the coherent evolution. Notice that the average fidelity 
quickly drops to small values specially when t g f c = 5. 

For a given a, the fidelity is less affected when the noise 
characteristic frequency is high (t g f c ^S> 1). In other 
words, the high-frequency noise generates less decoher- 
ence. This observation in our simulations can be justi- 
fied first of all in terms of our noise model specially by 
looking at the power spectral density given by Eq. (JoJ) . 
As the noise frequency increases (r c — > 0), the power 
spectral density function approaches zero for all u). This 
means that the noise disappears and the fidelity remains 
unaffected. 

However, an alternative way to analyze this result is 
using the theorem given in Ref. 7i. The time evolution 
of the system over each control time T = t c N is given by 



N 



Un(T) = [] exp[-*(H + e^H )r c ], 



(10) 



k=l 



where Ho and H are obtained from Hq and H, Eq. (fij) 
and Eq. ^ respectively, after replacing {Jij} by {JijJ- 
Expanding the product, we can write the result up to the 
first power of T as 



- N 

iHT ~ e 

k=l 



(fc) 



N 
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FIG. 1. (Color online) The average fidelity versus standard 
deviation a for t g f c = 1,2,5,10,20,40,100,200. The set of 
control fields used in each case has been obtained by maxi- 
mizing Eq. |5J, which leads to a fidelity about 99.83% with 
respect to the perfect system. 



But according to the weak law of large numbers 



lim 



N 



k=l 



> e = 0, 



and thus the first power of HqT as well as all the higher 
powers approach zero in probability, hence 



N 



lim (|| U N (T) - exp(-iHT) ||> e) = 0. 



(11) 



Therefore, the effect of noise disappears if the character- 
istic frequency t~ 1 is sufficiently high. 

The validity of Eq. (Ill can be also verified by de- 



creasing the number of realizations, that is, performing 
an average with smaller samples. In this case, the result- 
ing plots will be no longer smooth. Actually, the points in 
each plot are scattered around the average curve. How- 
ever, the smoothness of the plots corresponding to high 
frequencies are less affected. In another words, the plots 
for higher frequencies are less sensitive to the number of 
realizations. 

However, using the central limit theorem, r; = 
lim^v^^ X^feLi e ^ /V^V goes to a Gaussian distribution 
with the same mean and variance of e". Therefore, for 
finite but large N [7] 



(12) 



U N (T) ~ ex.p(-iHT)exp(-ir}H T/VN) 



which shows how the error enters in each control step 
T of the implementation and the errors will eventually 
reduce the gate fidelity. This result is not valid for in- 
termediate values of N, since the terms corresponding to 



the nonzero commutator brackets of Hq and H become 
also important. 

Fig. [T] also shows that for a fixed a the average fi- 
delity decreases when we decrease the noise frequency 
from tgf c = 200 until t g f c — 5 and the fidelity has the 
opposite behavior from t g f c = 5 until t g f c = 1. This 
behavior can be explained from the effect of the control 
pulses on an imperfect system. When the control-pulse 
frequency is close to the noise frequency, t^ 1 = T _1 , 
its effect helps to protect the system against decoher- 
ence, similar to dynamical decoupling [31] When we 
increase the noise frequency, the control-field effect be- 
comes less relevant, and for t g f c > 5 it plays almost no 
role at all. 

Finally, Fig. [T] shows that the curves apparently con- 
verge to the same average fidelity for sufficiently large 
(j. The limiting average fidelity is close to 1/d = 0.125, 
which corresponds to the one obtained from the full 
randomizing map e(p) = l/<i, where p is an arbitrary 
(pure) state and d the Hilbert space dimension. How- 
ever, this saturation value (0.125) can be obtained if the 
randomness in the gate parameters uniformly generates 
all elements in SU(8), and if we end up in a Lie sub- 
group of SU(8), the saturation value can be different (see 
Ref. [22]). What we observe after performing simulations 
with large values of a is that the saturation values are 
not exactly the same, but slightly different depending on 
the noise frequency. Therefore, fixing the control Hamil- 
tonian and just changing the coupling J does not seem 
to generate the whole SU(8) space. 

Now, we come back to justify why introducing one 
source of noise to J is sufficient to obtain essentially the 
same results of having three independent sources of noise 
affecting three couplings { J12, J13, </23j-- In the later case, 
the exponent in Eq. ( 10 1 is replaced by 



-i(H 



H )i 



where 



V fc 12 > fc 13 ' fc 23 I 



H n 



(Hi2,H 13 ,H 23 ), and 



By using ij — lim^^oo X^feLi e ^ fe V V^V instead of 77, a 
similar reasoning to the one used for one J can be ap- 
plied here to reach the same results as before for high 
noise frequencies. Moreover, we have repeated the simu- 
lations with three independent sources of noise and found 
essentially the same results for low noise frequencies as 
well. Actually, increasing the number of noise sources 
just leads to faster decay of average fidelity for all fre- 
quencies. Simulations with six different couplings (dif- 
ferent couplings in X and Y directions) and with six in- 
dependent sources of noise confirm the latter statement. 
In the new simulations, there is no universal saturation 
value. 
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C. Static Imperfections 



III A the 



As we have discussed at the end of Sec. 
case t c — > oo corresponds to a fixed noise in the sys- 
tem, that is, the noise does not change in time. Such a 
static noise can be associated with, for example, inacces- 
sibility of measuring exactly the system parameters such 
as the coupling constants. Static noises do not lead to 
any decoherence, but introduce error into the gate im- 
plementation. It is known that entanglement dynamics 
in similar systems remain almost unchanged under the 
influence of uniform static noise when it is smaller than 
10% (see Refs. 0E]). 

Analyzing the error for the static-imperfections case 
is mathematically equivalent to the case of dynamical 
noise with r^ 1 = ij 1 . Again, we adjust the control 
Hamiltonian for the perfect system and use it for a sys- 
tem with static imperfections. The coupling is given by 
Eq. Q with a fixed random variable e(t) — e for the 
whole evolution. Considering a uniform distribution for 
e, we generate a large sample of couplings with a given 
half-width 8 and then calculate the corresponding fidelity 
using Eq. After obtaining the sample average and re- 
peating the procedure for many values of 8, we find how 
the fidelity changes as a function of the half- width. More- 
over, we have decreased the noise level interval to be more 
focused on the region with unaffected entanglement. 

The solid curve in Fig. [2] depicts the average fidelity 
as a function of the half-width 8. The number of real- 
izations per each fixed 8 is 100,000. The diagram shows 
that the average fidelity decreases when the half-width 
8 increases. Here, we have focused on half-widths less 
than 0.5 because the experimental upper bound for such 
superconducting qubit chain is quite below that [5]. Ac- 
tually for disorders such as 8 < 10%, the average fidelity 
remains above 97.83%. Other curves in Fig. [2] will be 
discussed in the next section. 



IV. IMPROVING THE ROBUSTNESS 
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FIG. 2. (Color online) The average fidelity versus half- width 8 
(static imperfections) for three different set of control fields: 
solid line corresponds to the fields obtained by maximizing 
Eq. ((HI), dotted line corresponds to the fields obtained by solv- 
ing the new optimization problem ( 14 1 and dashed line shows 
the result after solving problem ( 16 1 with the global solution 
of problem 



([14j) as its initial guess. The corresponding fi- 
delities are about 99.83%, 99.47% and 99.83% for the perfect 
system, respectively. Inset: a zoom on the region < 8 < 0.1. 



the resulting optimal solution has the best fidelity at that 
specific coupling irrespective of the fidelity at points for 
which J =/= J. However, with noisy couplings a more 
robust set of control fields should have high fidelities also 
in the couplings that deviate from J — J. To this end, 
we modify the objective functional such that it includes 
a range of coupling values [J — 8 J, J +8 J]. We solve the 
new optimization problem 



arg 



J+SJ 



F(u, J)w{J) dJ 



(14) 



J-SJ 



In Sec. |III| we have considered the set of control fields 
optimized for a fixed value of coupling and then applied 
that set to a system with imperfections. Actually, we 
have investigated there the performance of that specific 
set of fields in the presence of noise. However, it is pos- 
sible to find another set of control fields that are more 
robust to noise. In this section, we use the optimal con- 
trol techniques to find a different set of control fields 
representing an improved performance in systems with 
imperfections. 

The set of control fields which is used in Sec. Illll is the 
solution u of the optimization problem 



arg < max F 



(13) 



where F is given by Eq. But, the objective functional 
F (u, J) uses a single coupling value J = J, and therefore 



where w(J) is an appropriate weight function enabling 
us to put different stress on the points in the interval 
brought into the process of optimization. 

The solid line in Fig. [3] depicts the fidelity in terms of 
J I J for the fields optimized for a single val ue of coupling 
J = J obtained from solving problem (13). Inspired by 
this plot, we set the weight function such that it has 
small values in the vicinity of J / J — 1 and large values 
in the outermost points of the interval. In fact, by simply 
letting 



w(J) 



0, |J_i|<^ 



1, <Si < \ j -1| <<5 2 , 



(15) 



we can find a new set of control fields leading to higher fi- 
delities in points different from J/ J. Solving the new op- 
timization problem with 200 random initial guesses and 



(> 



choosing the fields with highest fidelity, we can approach 
to the global solution in this case. The dotted line in 
Fig. [3] depicts the corresponding fidelity in terms of J/ J 
for such set of fields with Si = 0.05, S2 = 0.15, and 
SJ = 0.15J. It can be seen that the new optimal fields, 
comparing with the previous ones, have smaller fidelities 
in the vicinity of J/ J = 1 but larger in other points. 

However, we can improve the fidelities of this new set of 
control fields in the vicinity of J/ J = 1. We set w(J) = 1 
and add some new terms smoothening the fidelity and its 
derivative around the peak by resorting to 



arg < max 



J+SJ 



J-SJ 



F -a 



dF 
~dJ 



d 2 F 



dJ 2 



dJ\, (16) 



where a and /3 are some positive constants. Actually, 
the set of con trol fields just obtained by the optimization 
problem ( 14 1 can be improved in the vicinity | J— J| < 0.1 



when used as the starting point in problem ( 16 1 with 
a = 0.001 and /3 = 0.003. The dashed line in Fig. [3] 
shows the corresponding fidelity versus J/ J in this case. 

The performance of our new sets of control fields in 
the presence of noise can be analyzed in the same way of 



Sec. Ill In the case of static imperfections, the result for 



the fields optimized according to problem ( 14 1 with the 
weight function (151 is shown by dotted line in Fig. [5] 
The dashed line in the same plot shows the situation 
for the fields optimized solving problem (16). Specially, 



for 10% disorder (S = 0.1) the corresponding fidelities 
are 97.83% (problem fDj)), 98.47% (problem (JTi|), and 
98.67% (problem |l6| jshowing 0.64% improvement in 
the second case and 0.84% improvement in the last case. 

For the dynamical-imperfections case, when the noise 
frequency is high, the resulting plots for the three sets 
of control fields remain almost the same because the ef- 
fect of control pulses plays no role in high noise frequen- 
cies. However, for low noise frequencies t g f c = 5,2,1 (see 
Sec. Ill B I the new sets of fields have better performances. 



V. SUMMARY AND DISCUSSION 

In this paper we considered the effect of imperfections 
on a recently established Toffoli gate realized in circuit 
QED with quantum control methods. The total Hamilto- 
nian is bilinear and we studied the effect of imperfections 
on the system Hamiltonian by introducing noise in the 
interqubit couplings. We showed that in the case of dy- 
namical imperfections the average fidelity is less sensitive 
to noise for high characteristic frequencies. Actually the 
effect of noise completely disappears when the noise fre- 
quency is sufficiently high. For static noise we showed 



that the fidelity decreases by 2% when the system af- 
fected by a uniform noise with half-width less than 10%. 

We also obtained two new sets of control fields by mod- 
ifying the objective functional in the original optimiza- 
tion problem considered in Ref. [TB]. We showed that 
these new sets of control fields are more robust when af- 
fected by static noise. 
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FIG. 3. (Color online) Fidelity versus J / J for three different 
set of control fields: the solid line corresponds to the fields 
obtained by maximizing Eq. jEJ, the dotted line corresponds 
to the fields obtained by solving the new optimization prob- 
lem ( |14[ ), and the dashed line shows the result by solving 
problem (16 1 with the global solution of problem 1 14 1 as its 



initial guess. Inset: a zoom on the peak. 



We considered independent distributions in different 
time intervals for random variables in the noise model 
Eq. ^8j. It is possible to extend this method simulating 
other types of noise by considering various dependencies 
between the random variables in different times. Conse- 
quently, the average fidelity may be improved in a more 
effective way having those sorts of dependencies. The 
techniques used in this work may be applied to other 
multi-qubit gates as well. 
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